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ABSTRACT 

Aims. We aim to study excitation of the observed ~ 5-min oscillations in the solar corona by localized pulses that are 
launched in the photosphere. 

Methods. We solve the full set of nonlinear one-dimensional Euler equations numerically for the velocity pulse propa- 
gating in the solar atmosphere that is determined by the realistic temperature profile. 

Results. Numerical simulations show that an initial velocity pulse quickly steepens into a leading shock, while the non- 
linear wake in the chromosphere leads to the formation of consecutive pulses. The time interval between arrivals of two 
neighboring pulses to a detection point in the corona is approximately 5 min. Therefore, the consecutive pulses may 
result in the ~ 5-min oscillations that are observed in the solar corona. 

Conclusions. The ~ 5-min oscillations observed in the solar corona can be explained in terms of consecutive shocks that 
result from impulsive triggers launched within the solar photosphere by granulation and/or reconnection. 
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1. Introduction 

Propagating acoustic waves are frequently seen in the so- 
lar corona as periodic variations of spectral line intensity 
(D e Moortel e t al. 120001 ?IUU^ Marsh et al. Lin e t 

al. [20051 (20061 Srivastava et al. [20081 Wang et al. [2009|l . 
As these waves are often observed within the frequency 
range corresponding to the acoustic waves in the solar pho- 
tosphere/chromosphere, this logically leads to the idea of 
penetration of the photospheric acoustic oscillations into 
the corona. However, the photospheric 5-min oscillations 
are evanescent in the gravitationally stratified solar atmo- 
sphere as their frequency is lower than the cut-off frequency 
(Lamb [1^051 Roberts IMHl Musielak et al. 120051) . Bel & 
Leroy (|1977p suggested that the cut-off frequency of the 
magnetic field-free atmosphere is lower for waves propagat- 
ing obliquely to the vertical direction. De Pontieu et al. 
(|2005p proposed that p-modes may be channeled into the 
solar corona along inclined magnetic field lines as a result 
of the decrease of the acoustic cut-off frequency. Mcintosh 
& Jefferies |2006j) found the observational justification of 
the modification of the cut-off frequency by inclined mag- 
netic field. As the magnetic field of active region loops is 
predominantly vertical in the photosphere/chromosphere, 
it is unclear how p-modes may penetrate into the coronal 
regions. 

There are two different types of drivers in the highly 
dynamic solar photosphere: oscillatory (e.g. p-modes) and 
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impulsive (e.g. granulation and/or explosive events due to 
magnetic reconnection). Both types of drivers may be re- 
sponsible for the observed dynamical phenomena in upper 
atmospheric regions. 

As a result of the rapid decrease of mass density with 
height, finite-amplitude high-frequency photospheric oscil- 
lations can quickly grow in their amplitudes and steepen 
into shocks, which by energy dissipation can lead to the 
chromosphcric heating (Narain & Ulmschneider .19901 119961 
Carlsson & Stein 119971 Ruderman 2Q06}. Lower- frequency 
waves, those with ~5 min period, are not good candidates 
for the chromosphcric heating (Narain & Ulmschneider 
fTOOOl) . 

It was found by HoUweg (|1982p that a localized pulse 
that is launched initially sets up a nonlinear wake which 
results in a trail of consecutive shocks. Such shocks were 
called rebound shocks by HoUweg (1982). 

The time interval between consecutive shocks is close 
to the period of the nonlinear wake. In the linear case, 
the wake oscillates with the acoustic cut-off frequency 
of the stratified atmosphere. Nonlinearity modifies the 
wave period of the wake, with many features of spicules 
which exhibit the periodicity of about 5-min (Murawski & 
Zaqarashvili 120101) . Then, these quasi-periodic shocks may 
lead to the oscillatory dynamics of coronal plasma, which is 
observed as intensity oscillations in coronal spectral lines. 

The aim of this paper is to study the role of rebound 
shocks, which are formed by an impulsive perturbation, on 
the observed 5-min oscillations in the solar corona. Here 
we consider the simplest hydrodynamic case, which can be 
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developed to more realistic magnetohydrodynamic model 
in future studies. 

This paper is organized as follows. The basic equations 
and the atmospheric model are described in Sect. [5] The 
numerical model and results of numerical simulations of 
impulsive photospheric driver are discussed in Sect. [31 This 
paper is concluded by a summary of the main results in 
Sect, m 

2. Basic model 

2.1. Hydrodynamic equations 

Our model system is taken to be composed of a 
gravitationally-stratified solar atmosphere that is described 
by one-dimensional (ID) Euler equations: 



dg d{gV) 



dt 



dy 



0, 



(1) 
(2) 
(3) 



dV ^dV dp 

'^ + '^^ = -d-y-''' 

dt dy ^ ^'^ dy 

Here g denotes the mass density, y is a vertical compo- 
nent of the flow velocity, p — k^gT/m is the gas pressure, 
7 = 5/3 is the adiabatic index, g = 272 m s~^ is the gravi- 
tational acceleration, T is the temperature, m is the mean 
particle mass and is Boltzmann's constant. 



2.2. The equilibrium 

We assume that at the equilibrium the solar atmosphere 
is settled in a static {V — 0) environment in which the 
pressure gradient force is balanced by the gravity, that is 

-f-«^0. (4, 

With the use of the equation of state we obtain the equi- 
librium gas pressure and mass density as 



Po(y) = Poo exp 



dy 



Qoiv) 



po{y) 
gMy) 



(5) 



Here A{y) — k-^T [y) / {nig) is the pressure scale-height and 
Poo is the gas pressure at the reference level, chosen here at 
y^ = 10 Mm. 

We adopt a realistic temperature profile T{y) for the so- 
lar atmosphere (Vernazza et al. 1981). In this profile T at- 
tains a value of about 5700 K at the top of the photosphere 
which corresponds to y — 0.5 Mm. At higher altitudes T{y) 
falls off until it reaches its minimum of 4350 K at y ~ 0.95 
Mm. Higher up T{y) grows gradually with height up to the 
transition region which is located at y ~ 2.7 Mm. Here 
T(y) experiences a sudden growth up to the coronal value 
of 1.5 MK at ?/ = 10 Mm. Having specified T{y) we can 
obtain mass density and gas pressure profiles with the use 
of Eq. p. 

Equations lU-Q are to be solved numerically for an im- 
pulsive perturbation launched initially within the solar pho- 
tosphere. The numerical simulations for a harmonic driver 
were already reported elsewhere (e.g. Erdelyi et al. 120071 
Fedun et al. 120091 and references therein). Therefore, it is 
justifiable to limit our study to the case of impulsively gen- 
erated waves. 



2.3. Linear approximation 

Neglecting all nonlinear terms in Eqs. ([I])-® leads to the 
classical Klein-Gordon equation (Rae and Roberts 119821 
Roberts 2004) 



d^Q ^d^Q 



dt^ 



dy'^ 



nlQ = 0, 



where Q{t, y) = V{t, y) ^ goi^VM / Qo{y)cl{y) and 



(6) 



(7) 



In the case of an isothermal atmosphere, the Klein-Gordon 
equation implements a natural frequency of the stratified 
medium, 0,^ = Cs/2A, called the cut-off frequency for acous- 
tic waves. A physical interpretation of the cut-off frequency 
is that waves of lower (higher) frequencies than J7c are 
evanescent (propagating). According to the theory devel- 
oped for the linear Klein- Gordon equation the initially 
launched pulse results in a leading pulse which propagates 
with the sound speed. This pulse is followed by the wake, 
which oscillates with ilc and gradually decays as time pro- 
gresses fLamb IT908I Rae and Roberts 1982). However, sit- 
uation is changed when one considers nonlinear terms in 
Eqs. ((I])-©- Then, a large-amplitude initial pulse may lead 
to consecutive shocks due to nonlinearity (Hollweg 11982|) . 
In the nonlinear regime, the wake may oscillate with the 
period, which differs from the linear cut-off period. Then 
the interval between consecutive shocks may depend on the 
amplitude of the initial pulse (Murawski and Zaqarashvili 
I2010p . We solve the fully nonlinear equations for a realis- 
tic atmospheric model numerically and the results are pre- 
sented in the next section. 



3. Numerical simulations for impulsively generated 
waves 

Equations (HJ-® are solved with a use of the code FLASH 
(Dubey et al. 2009). We set the simulation region as —0.5 
Mm < y < 29.5 Mm and at the numerical boundaries we fix 
in time all plasma quantities to their equilibrium values. In 
our studies we use an adaptive mesh refinement (AMR) grid 
with a minimum (maximum) level of refinement blocks set 
to 3 (9). In the simulations, we excite waves by launching 
at t — s the initial velocity pulse of a Gaussian profile 



V{y, t = 0) = exp 



iy - yoY" 



(8) 



with Av the amplitude of the initial pulse, j/o its initial 
position and w its width. We set and hold fixed yo — 0.5 
Mm and w = 0.1 Mm, but allow to vary in different 
simulations. 

At the initial stage of the wave evolution the initial pulse 
splits into counter-propagating waves. As a result of the 
rapid decrease of the equilibrium mass density, the upward 
propagating pulse grows in its amplitude and quickly steep- 
ens into a shock. Figure [1] illustrates the spatial profile of 
V{y) aXt = 250 s for A^ ^ 1 km s~^. At this time, the pulse 
reaches the level y = 23.5 Mm, while the secondary shock 
begins to form from the nonlinear wake at y = 2 Mm. Later 
on, this secondary shock propagates upwards, and then it 
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Fig. 1. Velocity (in units of 1 km s~^) profile vs height y 
(in units of 1 Mm) a.t t = 250 s for = 1 km s~^. Here 
y = Mm {y ~ 2.7 Mm) corresponds to the base of the 
photosphere (the transition region). 



is followed by the next shock. This process repeats in time 
until the perturbation energy finally subsides to zero. 

Figure [5] displays the temporal dynamics of the velocity 
that is collected at y = 20 Mm for the initial pulse ampli- 
tude of Av = 1 km s^^ (top panel) and A^ — 2 km s^^ 
(bottom panel). In the top panel, the arrival of the lead- 
ing shock to the detection point occurs at i ~ 250 s and 
the second pulse reaches the detection point at i ~ 570 s 
i.e. after ~ 320 s. Thus, the nonlinear wake in the chromo- 
sphere excites quasi-periodic upward propagating pulses, 
which enter the corona and may cause quasi-periodic in- 
tensity oscillations which are observed in the imaging data. 
For Av = 1 km s~^, which is close to the typical granular 
velocity, the interval between arrivals of neighboring shocks 
is near 5-min (Fig. [21 top). This means that the nonlinear 
wake in the realistic atmosphere and for the initial pertur- 
bations, corresponding to the solar granular velocity, ex- 
hibits about a 5-min wave period, i.e. it is longer than in 
the linear isothermal case (~ 3-min in the chromosphere). 
Time-signature for A^ = 2 km s~^ shows that the leading 
shock arrives to the detection point at i ~ 190 s, while the 
second shock reaches this point at i ~ 630 s, i.e. after ~ 440 
s (Fig. [21 bottom panel). It is clearly seen that the interval 
between arrival times of two consecutive shocks is longer 
for a larger amplitude of the initial pulse. 

4. Discussion and summary 

Frequently observed ^ 5-min oscillations in the solar corona 
are often explained by leakage of photospheric p-modes 
along inclined magnetic field. Vertically propagating acous- 
tic waves with 5-min period are evanescent due to the strat- 
ification of the solar atmosphere as chromospheric acoustic 
cut-off period is ~ 3 min (Roberts 12004). But, acoustic- 
gravity waves have a smaller cut-off frequency when they 
propagate with the angle about the vertical. This may al- 
low the photospheric 5-min oscillations to channel along a 
stratified chromosphere and penetrate into the corona (De 
Pontieu et aL [20051 Erdelyi et al.|2007, Fedun et al.,20 0_9 ). 
In order to increase the cut-off period from 3-min up to 5- 
min, the propagation angle (or magnetic field inclination) 
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Fig. 2. Time-signatures of V (in units of km s^^) collected at 
y = 20 Mm for the case of = 1 km s~^ (top panel) and 
Av = 2 km s~^ (bottom panel). Time is expressed in units of 1 
s. 



should be 6* 50° {cos9 w 3/5). Therefore, the leakage of 
p-modes may take place only in particular regions of the 
solar atmosphere, where the magnetic field is significantly 
inclined in the chromosphere. 

In this paper, we suggest an alternative mechanism to 
explain the observed oscillations in the solar corona, which 
is based on the rebound shock model of Hollweg (|1982p . 
We numerically solved the full set of Euler equations for 
the realistic VAL-C temperature profile and for a Gaussian 
velocity pulse launched within the photosphere. We found 
that velocity pulses, originating from granules or magnetic 
reconnection in the lower regions, lead to different responses 
of the chromosphere/transition region than the periodic 
acoustic waves resulting from p-modes. It must be men- 
tioned, that the energy of granular motions is higher than 
the energy of p-modes, therefore the impulsively triggered 
waves should have more power than those triggered by a 
periodic driver. 

The numerical simulations show that as a result of the 
rapid decrease of the equilibrium mass density the initial ve- 
locity pulse quickly steepens into a shock. The shock prop- 
agates into the corona, while the nonlinear wake is formed 
in the chromosphere due to the atmospheric stratification. 
This nonlinear wake leads to consecutive shocks as was first 
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shown by Hollweg (jl982p . The interval between the arrival 
times of two consecutive shocks depends on the amplitude 
of the initial pulse; a stronger pulse leads to longer inter- 
vals. The initial pulse with the granular velocity of 1 km 
s""'^ leads to ~ 5-min intervals between consecutive shocks. 
Therefore, the quasi-periodic arrival of consecutive shocks 
in the solar corona may cause the intensity oscillations with 
a period close to the interval between the shocks. 

We implemented a simple ID analytical model in order 
to avoid the propagation of acoustic oscillations with the 
angle to the vertical. We showed that purely vertically prop- 
agating pulses may lead to quasi 5-min oscillations in the 
corona due to consecutive shocks. Therefore, it is not nec- 
essary to invoke the propagation along inclined magnetic 
field in order to explain the observed 5-min periodicity in 
the corona. 

One dimensional propagation for acoustic waves is jus- 
tified for purely vertical magnetic field. In this case, acous- 
tic waves are in fact slow magneto-acoustic waves for low 
plasma /? ~ Cj/u^ < 1 and fast magneto-acoustic waves 
for high plasma f3 > 1. Here va is the Alfven speed. In the 
solar photosphere /3 is larger than unity, but it rapidly de- 
creases due to the mass density fall off with height (and con- 
sequent increase of the Alfven speed). It becomes smaller 
than unity in the chromosphere (Gary I200ip and tends to 
unity somewhere between the photosphere and the chromo- 
sphere and this surface should be thinner than the width 
of the chromosphere. The linear fast and slow magneto- 
acoustic waves are coupled near the level of /3 ~ 1 when 
they propagate obliquely to the magnetic field (Rosenthal 
et al. 120021 Bogdan et al. I2003p . However, these waves re- 
main purely acoustic for parallel propagation unless the 
tube dispersive effects are taken into account. Therefore, 
a pulse propagating along the vertical magnetic field may 
not feel the (3^1 surface. On the other hand, the acoustic 
wake, which is formed behind the pulse and oscillates along 
the magnetic field, may lead to the non-linear energy trans- 
fer into Alfven waves near the /3 ~ 1 region (Zaqarashvili 
and Roberts 2006, Kuridze and Zaqarashvili 12008). This 
effect may be revealed as far as one includes the magnetic 
field into numerical simulations. Then, a part of oscillation 
energy may be transformed into transverse oscillations near 
this region, but most of energy will remain in longitudi- 
nal oscillations due to thinness of /3 ~ 1 region. Therefore, 
/3 ~ 1 region may not significantly affect the formation 
of rebound shocks in the chromosphere and consequently 
quasi-periodic acoustic oscillations in the lower corona. 
However, a two-dimensional consideration and inclusion of 
magnetic field are necessary for the complete understanding 
of the proposed scenario. The first step in this direction was 
already done by Murawski and Zaqarashvili (|2010p in the 
modeling of spicule formation, but they considered a simple 
temperature profile covering only the chromosphere-corona. 
The plasma /3 was considered less than unity along the 
whole simulation region, therefore the effects of the (3^1 
region were absent in these simulations. We intend to con- 
sider the rebound shock model in the case of magnetic field 
and realistic temperature profile in the future. 

An important consequence of the rebound shock model 
is that the interval between consecutive shocks depends on 
the initial amplitude of pulse (see Fig. 2). The smaller am- 
plitude pulses lead to shorter interval between consecutive 
shocks with lower limit of 3-min, which is the linear acoustic 
cut-off period of the solar atmosphere. The interval is longer 



for stronger initial pulses being ~ 5-min for 1 km s^^ in the 
photosphere. Then the rebound shock model predicts the 
longer period oscillations above the regions of strong gran- 
ular power. The granulation is suppressed in strong mag- 
netic field regions of the photosphere (e.g. sunspots), there- 
fore the initial pulses should have smaller amplitudes there, 
leading to the oscillation at near cut-off frequency. Indeed, 
the strong magnetic field regions (sunspots, magnetic net- 
work cores) show predominantly 3-min oscillations in the 
chromosphere. Note that some observations also show the 
3-min oscillations in the solar corona above sunspots (De 
Moortel et al. I2002p and other magnetic structures (Lin 
et al. I2005P . These oscillations can be excited as a conse- 
quence of consecutive shocks due to chromospheric 3-min 
oscillations. 

As a result, we expect that ^ 5-min oscillations can 
be detected above the regions where the granular energy 
is significant i.e. the quiet Sun and surroundings of active 
regions/magnetic network. This is consistent with observa- 
tions. On the other hand, the detection of 5-min oscillations 
in the quiet corona is not an easy task due to a relatively 
weak intensity of coronal lines. However, a careful analysis 
still can be performed. But, one should keep in mind that 
the dynamics of acoustic waves in the field-free regions and 
along the vertical magnetic field could be quite different 
as it was discussed above. An initial acoustic pulse may 
be spread horizontally in field-free regions, while almost 
the whole energy would be guided along the field lines in 
magnetic structures. Therefore, the amplitude of intensity 
oscillations could be smaller in the quiet corona than near 
active regions and chromospheric network cores. However, 
the divergence of magnetic field and increased thermal con- 
duction may significantly weaken the amplitude of coronal 
oscillations, which seem to be quite strong and non-linear 
on Figs. 1 and 2. Then, the strong slow wave pulses may be- 
come almost linear once they penetrate into the corona as it 
is seen by observations. Our numerical model then requires 
the inclusion of magnetic field and the thermal conduction 
(at least, in the coronal part of the atmosphere). This will 
be done in future studies. 

It should be noted that the granular velocities may take 
values between 0.5-2 km s~^ with peak on 1 km s~^. As the 
interval between consecutive shocks strongly depends on 
the amplitude of the initial pulse, then the resulted coronal 
oscillations may take values between 4-7 min with a peak on 
5-min. Indeed, the wavelet analysis of coronal line images 
obtained by Hinode/EIS show that the oscillation power of 
coronal oscillations is concentrated at the period in a range 
of 4-6 min (Wang et al. I2009P , which is fully consistent with 
our theory. 

Our simulations were performed for an isolated pulse in 
order to show clearly the effect of rebound shocks. However, 
the solar photosphere is very dynamic, hence the initial 
pulse probably is followed by other pulses. The subse- 
quent pulse coming from the photosphere may also trigger 
the consecutive shocks in the chromosphere. The interac- 
tion between rebound shocks formed by different photo- 
spheric pulses may set up complex dynamics in the chro- 
mospheric plasma with immediate influence on the lower 
corona. Therefore, the coronal oscillations may have broad 
spectrum as we already discussed in the previous para- 
graph. On the other hand, if the mean interval between 
subsequent initial pulses is close to the mean interval be- 
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tween consecutive shocks, then a resonance may occur in 
the chromosphere. This could be subject of future study. 

The excitation of coronal oscillations due to leakage of 
p-modcs may occur only along significantly inclined mag- 
netic field in the chromosphere (in order to reduce the cut- 
off frequency). The magnetic field lines, which are signifi- 
cantly inclined from the vertical in the chromosphere, may 
not reach the corona at all. Therefore, the p-mode leak- 
age may have problems in real geometry of active region 
magnetic field. On the contrary, the rebound shock mecha- 
nism may work in any geometry of magnetic field including 
the purely vertical field lines. Therefore, the excitation of 
5-min oscillations in the solar corona by photospheric im- 
pulsive drivers has larger area of application than that of by 
p-modes. We believe that the future sophisticated models 
may shed light on the excitation of coronal acoustic waves. 

Our conclusions are: 

(a) a velocity pulse that is initially launched at the photo- 
spheric level (due to granules or reconncction) quickly 
steepens into a shock and can penetrate into the corona, 
while a nonlinear wake that is formed behind this shock 
leads to consecutive shocks in the chromosphere; 

(b) for the initial photospheric pulse amplitude of 1 km 
the time interval between two consecutive shocks is ~ 5- 
min; the consecutive shocks propagate upwards and may 
cause the observed 5-min intensity oscillations in the 
solar corona; 

(c) the final conclusion is that the observed ^ 5-min oscil- 
lations in the solar corona could be caused by impulsive 
photospheric perturbations (convection, reconnection) 
not necessarily by p-modes. 
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